Solutions of renormalization group flow equations with full momentum dependence 
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We demonstrate the power of a recently-proposed approximation scheme for the non-perturbative 
renormalization group that gives access to correlation functions over their full momentum range. 
We solve numerically the leading-order flow equations obtained within this scheme, and compute 
the two-point functions of the O(N) theories at criticality, in two and three dimensions. Excellent 
results are obtained for both universal and non-universal quantities at modest numerical cost. 
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The renormalization group, in its non-perturbative 
version [l|, Q (also referred to as the exact renormal- 
ization group), provides a general formalism giving ac- 
cess, for arbitrary coupling strength, to a whole set 
of physically important quantities, universal as well 
as non- universal [H, thermodynamic functions and 
momentum-dependent correlation functions, etc. How- 
ever, most studies within this framework involve ap- 
proximations that restrict their scope to the calculation 
of thermodynamical quantities, or correlation functions 
with vanishing external momenta. In order to access 
the full momentum dependence, Blaizot, Mendez-Galain, 
and Wschebor (BMW) have recently introduced an ap- 
proximation scheme which overcomes this limitation[5]. 
In principle, this scheme allows to compute, in all di- 
mensions, at and away from criticality, both univer- 
sal and non-universal quantities, as well as momentum- 
dependent properties from p = up to the ultra-violet 
cut-off A (inverse lattice spacing). 

In this Letter, we present the first complete implemen- 
tation of the leading order approximation of the BMW 
scheme, and demonstrate its power by using O(N) mod- 
els as a testbed. We compute the entire momentum de- 
pendence of the two-point functions in two and three di- 
mensions and obtain excellent results for both universal 
and non- universal quantities. 

We start by a brief outline of the formalism. In or- 
der to simplify the presentation, we shall write only the 
equations corresponding to the case of a scalar field the- 
ory with quartic coupling, i.e., restrict the presentation 
to the case N — 1 (corresponding to the Ising model). 
The strategy of the renormalization group is to build a 
family of theories indexed by a momentum scale param- 
eter k, such that fluctuations are smoothly taken into 
account as k is lowered from the microscopic scale A 
down to 0. In practice, this is achieved by adding to 
the original Euclidean action S a mass-like term of the 
form ASkM = h S q Rk(q 2 MqM-q)- The cut-off func- 
tion R k (q 2 ) is chosen so that Rk(q 2 ) ~ k 2 for q < k, 
which effectively suppresses the modes ip(q < k), and so 



that it vanishes for q > k, leaving the modes ip(q > k) 
unaffected. One then defines a scale-dependent partition 
function 



Z k [J] 



Vip e 



-S[<p]-AS h [<p]+JJ<p 



(1) 



and a scale-dependent effective action Tk[4>] through a 
(slightly modified) Legendre transform 



\ogZ k [J] = jj<t>~\J Rk(q 2 )M- q > 



(2) 



with <j) = 6 In Zk/5J. The variation of the effective action 
Tk[<j>] as k varies is governed by Wetterich's equation^]]: 



d k T k = ~J d k R k (q 2 )G k [q, 



(3) 



-.(2)r 



where G k [q,<f>] = (T k >[q, <p] + R k {q 2 ))-\ and lj'[q,(p\ is 
the second functional derivative of Tk[<f>] w.r.t. <p. The 
initial conditions of the flow equation ^ correspond to 
the microscopic scale k = A where all fluctuations are 
frozen by the AS k term, so that Fj. = a[0] = S [</>]. The 
effective action of the original theory is obtained as the 
solution of © for k — > where R k (q 2 ) vanishes. Differ- 
entiating Eq. © m times with respect to <f> yields the 
flow equation for the vertex function T k [q\, . . . , q m ; <f>] . 
Thus for instance, the flow equation for I^ 2 ) reads: 



W(p) = J q d k R k (q 2 )G 2 (q) [T k d >(p-p-q,q)x 
G k (p+q)T k 3 \-p,p+q-q)~±r ( k 4:) (p,-p,q,-q) 



(4) 



(Here we assumed the field <fi to be uniform, and omitted 
the (f) dependence to alleviate the notation.) Note that 



the flow equation for T 



(m) 



.Qi, 



.q n 



involves T 



m+l 



and p^ 1 ™ -1-2 ) ; leading to an infinite hierarchy. 

The flow equation ([3]), and the equivalent flow equa- 
tions for the vertex functions, are exact, but their solu- 
tion requires, in general, approximations. It is precisely 
one of the virtues of this formulation of field theory to 



2 



suggest approximation schemes that are not easily de- 
rived in other, more conventional approaches. In par- 
ticular, one can develop approximation schemes for the 
effective action itself, that is, which apply to the entire 
set of correlation functions. The BMW approximation [Hj 
is such a scheme. It relies on two observations. First, 
the presence of the cut-off function Rk{q 2 ) insures the 
smoothness of the 's and limits the internal momen- 
tum q in equations such as Eq. |4]) to q < k. In line with 
this observation, one neglects the q-dependence of the 
vertex functions in the r.h.s. of the flow equations (e.g. 
in I^ 3 ) and in Eq. (01), while keeping the full depen- 
dence on the external momenta pi. The second observa- 

(m+l)j 

which enables one to close the hi- 



tion is that, for uniform fields, r^'" T - L "' (p 1: . . . } p m , 0, <j>) 



,Pm,' 

erarchy of equations. 

At the leading order of the BMW scheme one keeps 
the non trivial momentum dependence of the two-point 
function and implements the approximations above on 
Eq. |4|), which becomes: 



fcdfcifW) 

with 



J: 



r 



(2) 



d 2 A 2) (5) 



J n (pA)= / kdkRkiq^G^iq^^Gkip+q,^) . (6) 

Jg 

The approximation can be systematically improved: The 
order m consists in keeping the full momentum depen- 



dence of r 

r (m+2) 



(2) 
k ' ' 



r[ an d truncating that of V 



(m+1) 



and 



along the same lines as those leading to Eq. 
corresponding to m = 2. 

The zeroth order approximation is the so-called local 
potential approximation (LPA) where vertex functions 
are obtained as derivatives of the effective potential 14 
(equal, to within a volume factor, to I\. evaluated for a 



uniform 



^(m) 



(Pi 



V^ m \(j)), except for 



T^'(p,<jt) ^=p 2 +V^ ] (</>). The LPA has been widely used 
with reasonable success 0, 0, @, ■ It can be improved 
through a systematic expansion in gradients of the fields, 
usually referred to as the derivative expansion (DE) 0,0. 
However, in contrast to the BMW scheme, the DE, at any 
finite order, does not give access to correlation functions 
with non-vanishing external momenta (or with external 
momenta larger than the smallest mass). 

We now turn to the main purpose of the present Let- 
ter, which is to show that the nonlinear integro-partial- 
differential equation ((5l) can be studied as is, without 
further approximation J9|. Note that the earlier studies 
of Eq. ([5]) presented in [ll| involve additional approxima- 
tions which are linked to a specific cut-off function, and 
which become too crude below three dimensions. 

In order to treat efficiently the low (including zero) mo- 
mentum sector, we work with dimensionless and renor- 
malized quantities. Thus, we measure all momenta in 



units of k: p — p/k. We also rescale p = ^(f> 2 accord- 
ing to p = k 2 - d Z k K^p (K d = (2ir)- d S d /d, S d be- 
ing the volume of the unit sphere), and set 1^ (p,p) = 

k~ 2 ZZ^T^ (p, p) . The running anomalous dimension i] k 
is defined by kd k Z k = —j] k Z k , so that at a fixed point 
Zk ~ k~ v , r\ being the anomalous dimension of the field 
at the fixed point. The absolute normalization of Z k is 



(2) 



1. 



fixed by choosing a point (po, po) where dpT 

~ (2) 

Then, the flow equation of F^ {p, p) follows trivially from 
Eq.([5]). For numerical reasons, we actually solve two 



equations: one for (p, p) = p~ 2 [T ^\p, p) — T^ 1 (0, p)) — 
1 and one for the derivative of the dimensionless ef- 
fective potential Wk(p) = ^ 1 fc _2 p Vfe(p). Note that 
f^ 2) (0,p) = W k {p) + 2pWl{p). (Here and below, primes 
denote derivative w.r.t. p.) These two equations read 
(dropping the k index to simplify the notation): 

d t Y = 7 lk (l + Y)+pd p Y -(2-d-r lk )pY' 



-2pp 



(p 2 Y'+x k yj 3 (p,p)-xij 3 (o,p) 



-M0,p)(Y'/2 + pY") 



(7) 



d t W 



{ Vk -2)W+ (d~2 + Vk )pW' + -J[(0,p). (8) 



Here dt = kd k , r\ k is obtained by setting Y k [po,po] — 
in Eq.©, J„(p,p) = K^Z^k^-^J^p) and 
\ k (p)=3Wl(p) + 2pW%(p). 

In practice, we use a fixed, regular, (p, p) grid and 
restrict the range of the cut-off function by setting 
R k {q > 4) = 0. When computing the double integrals 
J 3 {p,p), we need to evaluate Y for momenta p + q be- 
yond pmax, the maximal value on the grid. In such cases, 
we set Y(p) ~ y(j) max ), an approximation checked to 
be excellent for p max > 5. To access the full momen- 

(2) 

turn dependence, we also calculate T k (p, p) at a set 
of fixed, freely chosen, external p values. For a given 
such p, p/k is within the grid at the beginning of the 
flow. This is no longer so when k < p/pma,^, then, we 
switch to the dimensionful version of ([7]), and also set 
JsiPt P) = G(p, p)J2(0, p), an excellent approximation 
when p > fcp max . 

We found that the simplest time-stepping (explicit Eu- 
ler), a finite-difference evaluation of derivatives on a reg- 
ular (p, p) grid, and the use of Simpson's rule to calcu- 
late integrals, are sufficient to produce stable and fast- 
converging results. For all the quantities calculated, the 
convergence to three significant digits is reached with a 
(p, p) grid of 50 x 60 points; with such a grid, a typical 
run takes a few minutes on a current personal computer. 

Physical quantities exhibit a small dependence on the 
shape of R k (q 2 ) and on the point (po,Po) where rj k is 
computed. Since in the absence of any approximation, 
they would be strictly independent of the cut-off func- 
tion and of the choice of the renormalization point, a 
study of this spurious dependence provides an indica- 
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FIG. 1: (Color online) Typical results at criticality and k = 
for TV = 2 in d — 3 (p and u are measured in units of 
A, u = 37r 2 10 -5 , r c ~ 6.802854731032857, a = 2.25). (a) 
dimensionful function Y(p,0) + 1 = T (2) (p, 0)/p 2 ; (b) same 
data as (a). Top panel: T (2) (p, 0)/p 2 ~ p"^ with 77 ~ 0.041 
for p — > (red dashed line). Bottom panel: expected UV 
scaling for the self-energy; the dashed line shows the exact 
two-loop result for the slope. 



tion of the quality of the present approximation. To 
this end, we use the family of cut-off functions Rk(q 2 ) = 
aZk q 2 / (exp(q 2 / k 2 ) — 1) and vary systematically the pa- 
rameters a, poi an( i Pq- I n au cases studied, we find the 
dependence on po and po to be much smaller than that 
on a, so that only the latter needs to be considered. As a 
function of a, physical quantities typically exhibit a sin- 
gle extremum a*, located near a = 2, which moreover, 
always points towards the best numerical estimates. Fol- 
lowing the principle of minimal sensitivity (PMS) [l2T ] , we 
regard these extremum values, being locally independent 
of a, as our best values. 

We now turn to the discussion of results obtained, at 
criticality, first in dimension d = 3 and for various val- 
ues of N. The initial condition of the flow is taken to 
be r^' (p, p) = p 2 + r + up, where the bare coupling 
u sets a scale, [u] = [p], independent of the cut-off A. 

Keeping u fixed, the critical value r = r c is found by 

(2) 

dichotomy. All expected features of T k (p) at critical- 
ity are observed, as illustrated in Fig. [TJ Fig. QJi shows 
the typical shape of Y(p, 0). In the infrared (IR) regime 
k < p « u, 4 2) (p, 0) = P 2 (Y(p, 0) + 1) - p 2 -\ (Fig. rjb, 
upper panel) . This IR behavior of 1^ (p) can be used 
to extract the value of 77; the value thus obtained is in 
excellent agreement with that deduced directly from the 
renormalization condition. The ultraviolet (UV) regime 
k, u <C p "C A exists if u is sufficiently small; this regime 
can be studied perturbatively and one finds that, in lead- 
ing order, p 2 F(p, 0) ~ it 2 log(p/w). The present ap- 
proximation reproduces this logarithmic behaviour with, 
however, a prefactor 8% larger than the two-loop result 
(Fig. [TJ), lower panel). Note that the complete two- loop 



behavior can be recovered by a simple improvement of 
the BMW scheme [ijj. 

A quantity particularly sensitive to the UV-IR 
crossover region is the shift, due to interactions, of the 
critical temperature of the dilute Bose gas In the 

limit of small coupling, the shift is proportional to u. The 
proportionality coefficient is given (with the normaliza- 
tion used in Ref. [13j |) by the non- universal quantity 



d 3 p 



1 



1 



r( 2 )(p) p 2 



(9) 



in the limit u —* 0. Note that the integrand is peaked 
at values of p ~ u, and is significant for momenta typi- 
cally in the range [iVtt/100, 10Nu]. Initially introduced 
for N — 2, corresponding to Bose-Einstein condensation, 
c is often used as a sensitive benchmark of various ap- 
proximations, as it tests the 2-point function over a wide 
range of momenta. It has been computed, for several 
values of N, on the lattice and with high-order (six-loop) 
perturbation theory (see Table HJ . 

Table Q] contains our results for c and the critical ex- 
ponents ij, v and ui, together with some of the best es- 
timates available in the literature. Our numbers are all 
given for the PMS values a* of the cut-off parameter, 
and the digits quoted remain stable when a varies in the 
range [a* — |, a* + |], The quality of these numbers 
is obvious: For all N values where six-loop resummed 
calculations exist, our results for c are within the error 
bars (and comparable to those obtained from an approx- 
imation specifically designed for this quantity [l6|); the 
results for v agree with previous estimates to within less 
than a percent, for all N; as for the values of 77 and ui, they 
are typically at the same distance from the Monte-Carlo 
and temperature series estimates as the results from re- 
summed perturbative calculations. For N = 100, we find 
c = 2.36, 77 = 0.0023, and v = 0.990, which compare 
well to the exact large N value c ~ 2.33 [13| and to the 
values 77 ^0.0027 and v = 0.989 obtained in the 1/N ex- 



pansion 15j. Our numbers also compare favorably with 
those obtained at order d 2 in the DE scheme 0. 

The two-dimensional case, for which exact results ex- 
ist, provides an even more stringent test of the BMW 
scheme. We focus here on the Ising model N — 1 which 
exhibits a standard critical behavior in d = 2, and the 
corresponding critical exponents. (The coefficient c is not 
defined in d = 2.) The perturbative method that works 
well in d — 3 fails here: for instance, the fixed-dimension 
expansion that provides the best results in d = 3 yields, 
in d = 2 and at five loops, 77 = 0.145(14) j2?| in contra- 
diction with the exact value rj = 4 [29( . We find instead 
rj = 0.254, v = 1.00, and u = 1.28 in excellent agreement 
with the exact values 77 = 4, v = 1 and the conjectured 
value to = § [H. 

To summarize, our results show that the single equa- 
tion §5§ (and its generalization to O(N) models) is suf- 
ficient to obtain the momentum dependence of the two- 
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TABLE I: Coefficient c and critical exponents of the O(N) models for d = 3. 



N BMW Resummed perturbative expansions 

rj v u c r\ v u c Ref.° 



Monte-Carlo and high-temperature series 
77 v lo c Ref. a 



0.034 0.589 0.83 

f 0.039 0.632 0.78 1.15 

2 0.041 0.674 0.75 1.37 

3 0.040 0.715 0.73 1.50 

4 0.038 0.754 0.72 1.63 
10 0.022 0.889 0.80 



0.0284(25) 0.5882(11) 0.812(16) [17] 

0.0335(25) 0.6304(13) 0.799(11) 1.07(10) [17] [14] 

0.0354(25) 0.6703(15) 0.789(11) 1.27(10) [17] [14] 

0.0355(25) 0.7073(35) 0.782(13) 1.43(11) [17] [14] 

0.035(4) 0.741(6) 0.774(20) 1.54(11) [17] [14] 

0.024 0.859 [26] 



0.030(3) 0.5872(5) 0.88 [18] [19] 

0.0368(2) 0.6302(1) 0.821(5) 1.09(9) [20] [21] 

0.0381(2) 0.6717(1) 0.785(20) 1.32(2) [22J[23J 

0.0375(5) 0.7112(5) 0.773 [24, 25] 

0.0365(10) 0.749(2) 0.765 1.6(1) [25] [21] 



"The first reference is for the critical exponents, the second for c. 



point function with excellent accuracy, in all momentum 
regimes, for all N, and in any dimension. All this is ob- 
tained at a modest numerical cost using simple numerical 
techniques. The study presented here is only the leading- 
order of a systematic approximation scheme. A study of 
the higher orders would be necessary in order to quan- 
tify the accuracy that has been reached. However, the 
robustness of our results can already be gauged from the 
weak residual dependence on the cut-off function. 

We focused here on critical theories since numerous 
and accurate results exist for the critical regime, allowing 
for detailed and systematic checks, but it is clear that the 
method can be also used to deal with generically simpler 
situations. For instance, one could calculate the structure 
factor as a function of the momentum and the correlation 
length, which is of experimental interest. The effect of 

an external magnetic field could also be investigated by 

(2) 

taking advantage of the built-in field dependence of T k . 
A detailed investigation of the d — 2, A" > 1 cases, is also 
at hand. Finally, this approach is not limited to O(N) 
theories. It can also be applied to disordered, nonequi- 
librium, or quantum systems, expanding from existing 
studies within the DE scheme of, e.g. absorbing phase 
transitions 0, 0] , or random-field models [3(| • 
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